************************************************************
* File Name: PartisanWorldviews_Replication.do
* Created by: Rebecca Kreitzer
* Created on: 1/24/22
* Purpose:	JOP Replication File
* Last modified on: 1/24/22
* Last modified by: Rebecca J. Kreitzer
************************************************************

clear

********************************************
* Figure 1. Social Constructions of Target Populations Mapping
* this is descriptive data visualized in the program Tableau
********************************************


********************************************
* Figure 2. Histogram: Distribution of "Consensus"
*(Standard Deviations of Deservingness Scores)
********************************************

use "Exported Data Files/sd_allgroups_sdcat", clear

hist des_all_sd_sdcat2, bin(8) ///
	graphregion(color(white)) xtitle("SD of SD for Deserving (All Respondents)") ///
	color(gs6) fcolor(gs5) 
	
	graph save Figures/hist_des_all, replace
	
********************************************
* Figure 3. High Consensus Groups by Partisanship
********************************************
use "Exported Data Files/means_sdcats_consensus_merged.dta", clear

gen mean_des_dem_flip = 100-mean_des_dem
gen mean_des_rep_flip = 100-mean_des_rep

twoway scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "asian", ///
		msize(small) msymbol(circle) ylabel(none) xlabel(none)  leg(label(1 "Asian (D)") pos(3) col(2), , size(vsmall)) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "asian", ///
		msize(small) msymbol(circle) ylabel(none) xlabel(none) leg(label(2 "Asian (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "workingclass", ///
		msize(small) msymbol(circle_hollow) ylabel(none) xlabel(none)  leg(label(3 "Working Class (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "workingclass", ///
		msize(small) msymbol(circle_hollow) ylabel(none) xlabel(none) leg(label(4 "Working Class (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "latino", ///
		msize(small) msymbol(smcircle) ylabel(none) xlabel(none)  leg(label(5 "Latinos (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "latino", ///
		msize(small) msymbol(smcircle) ylabel(none) xlabel(none) leg(label(6 "Latinos (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "middleclass", ///
		msize(small) msymbol(smcircle_hollow) ylabel(none) xlabel(none)  leg(label(7 "Middle Class (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "middleclass", ///
		msize(small) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(label(8 "Middle Class (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "bluecollar", ///
		msize(small) msymbol(square) ylabel(none) xlabel(none)  leg(label(9 "Blue Collar (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "bluecollar", ///
		msize(small) msymbol(square) ylabel(none) xlabel(none) leg(label(10 "Blue Collar (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "singledad", ///
		msize(small) msymbol(square_hollow) ylabel(none) xlabel(none)  leg(label(11 "Single Dads (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "singledad", ///
		msize(small) msymbol(square_hollow) ylabel(none) xlabel(none) leg(label(12 "Single Dads (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "africanamericans", ///
		msize(small) msymbol(smsquare) ylabel(none) xlabel(none)  leg(label(13 "African Americans (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "africanamericans", ///
		msize(small) msymbol(smsquare) ylabel(none) xlabel(none) leg(label(14 "African Americans (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "singlemom", ///
		msize(small) msymbol(smsquare_hollow) ylabel(none) xlabel(none)  leg(label(15 "Single Moms (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "singlemom", ///
		msize(small) msymbol(smsquare_hollow) ylabel(none) xlabel(none) leg(label(16 "Single Moms (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "mothers", ///
		msize(small) msymbol(triangle) ylabel(none) xlabel(none)  leg(label(17 "Mothers (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "mothers", ///
		msize(small) msymbol(triangle) ylabel(none) xlabel(none) leg(label(18 "Mothers (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "teachers", ///
		msize(small) msymbol(triangle_hollow) ylabel(none) xlabel(none)  leg(label(19 "Teachers (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "teachers", ///
		msize(small) msymbol(triangle_hollow) ylabel(none) xlabel(none) leg(label(20 "Teachers (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "youngblackmen", ///
		msize(small) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(label(21 "Young Black Men (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "youngblackmen", ///
		msize(small) msymbol(smtriangle) ylabel(none) xlabel(none) leg(label(22 "Young Black Men (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "jews", ///
		msize(small) msymbol(smtriangle_hollow) ylabel(none) xlabel(none)  leg(label(23 "Jews (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "jews", ///
		msize(small) msymbol(smtriangle_hollow) ylabel(none) xlabel(none) leg(label(24 "Jews (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "homeowners", ///
		msize(small) msymbol(diamond) ylabel(none) xlabel(none)  leg(label(25 "Homeowners (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "homeowners", ///
		msize(small) msymbol(diamond) ylabel(none) xlabel(none) leg(label(26 "Homeowners (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "millennials", ///
		msize(small) msymbol(diamond_hollow) ylabel(none) xlabel(none)  leg(label(27 "Millennials (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "millennials", ///
		msize(small) msymbol(diamond_hollow) ylabel(none) xlabel(none) leg(label(28 "Minnennials (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "native", ///
		msize(small) msymbol(smdiamond) ylabel(none) xlabel(none)  leg(label(29 "Native Americans (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "native", ///
		msize(small) msymbol(smdiamond) ylabel(none) xlabel(none) leg(label(30 "Native Americans (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_high == 1 & group == "retirees", ///
		msize(small) msymbol(smdiamond_hollow) ylabel(none) xlabel(none)  leg(label(31 "Retirees (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_high == 1 & group == "retirees", ///
		msize(small) msymbol(smdiamond_hollow) ylabel(none) xlabel(none) leg(label(32 "Retirees (R)")) mcolor(gs10) ///	
	graphregion(color(white)) ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(50, lwidth(vthin) lcolor(gs15)) ///
	ysc(r(0 100)) xsc(r(0 100)) ///
	xtitle(Deserving) ytitle(Power) ///
	aspectratio(1) ///
	subtitle("High Consensus") ///
	text(-5 10 "dependent", size(vsmall)) ///
	text(-5 90 "deviant", size(vsmall)) ///
	text(103 10 "advantaged", size(vsmall)) ///
	text(103 90 "contender", size(vsmall)) 
	
	saving(Figures/means_partisan_highcons, replace)


********************************************
* Figure 4. Low Consensus Groups by Partisanship
********************************************

	
	twoway scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "nra", ///
		msize(small) msymbol(circle) ylabel(none) xlabel(none)  leg(label(1 "NRA (D)") pos(3) col(2), , size(vsmall)) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "nra", ///
		msize(small) msymbol(circle) ylabel(none) xlabel(none) leg(label(2 "NRA (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "ice", ///
		msize(small) msymbol(circle_hollow) ylabel(none) xlabel(none)  leg(label(3 "ICE (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "ice", ///
		msize(small) msymbol(circle_hollow) ylabel(none) xlabel(none) leg(label(4 "ICE (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "whitenat", ///
		msize(small) msymbol(smcircle) ylabel(none) xlabel(none)  leg(label(5 "White Nationalists (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "whitenat", ///
		msize(small) msymbol(smcircle) ylabel(none) xlabel(none) leg(label(6 "White Nationalists (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "gunmanuf", ///
		msize(small) msymbol(smcircle_hollow) ylabel(none) xlabel(none)  leg(label(7 "Gun Manufacturers (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "gunmanuf", ///
		msize(small) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(label(8 "Gun Manufacturers (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "bigcorps", ///
		msize(small) msymbol(square) ylabel(none) xlabel(none)  leg(label(9 "Big Corporations (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "bigcorps", ///
		msize(small) msymbol(square) ylabel(none) xlabel(none) leg(label(10 "Big Corporations (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "bigbanks", ///
		msize(small) msymbol(square_hollow) ylabel(none) xlabel(none)  leg(label(11 "Big Banks (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "bigbanks", ///
		msize(small) msymbol(square_hollow) ylabel(none) xlabel(none) leg(label(12 "Big Banks (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "unauthimm", ///
		msize(small) msymbol(smsquare) ylabel(none) xlabel(none)  leg(label(13 "Unauth. Immigrants (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "unauthimm", ///
		msize(small) msymbol(smsquare) ylabel(none) xlabel(none) leg(label(14 "Unauth. Immigrants (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "richpeople", ///
		msize(small) msymbol(smsquare_hollow) ylabel(none) xlabel(none)  leg(label(15 "Rich People (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "richpeople", ///
		msize(small) msymbol(smsquare_hollow) ylabel(none) xlabel(none) leg(label(16 "Rich People (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "illegalaliens", ///
		msize(small) msymbol(triangle) ylabel(none) xlabel(none)  leg(label(17 "Illegal Aliens (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "illegalaliens", ///
		msize(small) msymbol(triangle) ylabel(none) xlabel(none) leg(label(18 "Illegal Aliens (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "police", ///
		msize(small) msymbol(triangle_hollow) ylabel(none) xlabel(none) leg(label(19 "Police (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "police", ///
		msize(small) msymbol(triangle_hollow) ylabel(none) xlabel(none) leg(label(20 "Police (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "military", ///
		msize(small) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(label(21 "Military (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "military", ///
		msize(small) msymbol(smtriangle) ylabel(none) xlabel(none) leg(label(22 "Military (R)")) mcolor(gs10) || ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "congress", ///
		msize(small) msymbol(smtriangle_hollow) ylabel(none) xlabel(none)  leg(label(23 "Congress (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "congress", ///
		msize(small) msymbol(smtriangle_hollow) ylabel(none) xlabel(none) leg(label(24 "Congress (R)")) mcolor(gs10) || ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "insuranceco", ///
		msize(small) msymbol(diamond) ylabel(none) xlabel(none)  leg(label(25 "Insurance Companies (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "insuranceco", ///
		msize(small) msymbol(diamond) ylabel(none) xlabel(none) leg(label(26 "Insurance Companies (R)")) mcolor(gs10)|| ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "opioiduser", ///
		msize(small) msymbol(diamond_hollow) ylabel(none) xlabel(none)  leg(label(27 "Opioid User (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "opioiduser", ///
		msize(small) msymbol(diamond_hollow) ylabel(none) xlabel(none) leg(label(228 "Opioid User (R)")) mcolor(gs10) 	|| ///
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "welfarecheats", ///
		msize(small) msymbol(smdiamond) ylabel(none) xlabel(none)  leg(label(29 "Welfare Cheats (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "welfarecheats", ///
		msize(small) msymbol(smdiamond) ylabel(none) xlabel(none) leg(label(30 "Welfare Cheats (R)")) mcolor(gs10) 	|| ///	
	scatter mean_pow_dem mean_des_dem_flip if consensus_all_des_low == 1 & group == "gunowners", ///
		msize(small) msymbol(smdiamond_hollow) ylabel(none) xlabel(none)  leg(label(31 "Gun Owners (D)")) mcolor(gs2) || ///
	scatter mean_pow_rep mean_des_rep_flip if consensus_all_des_low == 1 & group == "gunowners", ///
		msize(small) msymbol(smdiamond_hollow) ylabel(none) xlabel(none) leg(label(32 "Gun Owners (R)")) mcolor(gs10) ///		
	graphregion(color(white)) ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(50, lwidth(vthin) lcolor(gs15)) ///
	ysc(r(0 100)) xsc(r(0 100)) ///
	xtitle(Deserving) ytitle(Power) ///
	aspectratio(1) ///
	subtitle("Low Consensus") ///
	text(-5 10 "dependent", size(vsmall)) ///
	text(-5 90 "deviant", size(vsmall)) ///
	text(103 10 "advantaged", size(vsmall)) ///
	text(103 90 "contender", size(vsmall)) 
	

	saving(Figures/means_partisan_lowcons, replace)

********************************************
* Figure 5. Partisan Assessment of Low Consensus Groups
********************************************
use PartisanWorldviews_Data_Final, clear

twoway scatter pow_nra des_nra_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_nra des_nra_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("NRA (Pow.)") xtitle("NRA (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/nra_si, replace	
 
twoway scatter pow_ice des_ice_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_ice des_ice_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("ICE (Po.)") xtitle("ICE (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/ice_si, replace	

twoway scatter pow_whitenat des_whitenat_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_whitenat des_whitenat_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("White Nation. (Pow.)") xtitle("White Nation. (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/whitenat_si, replace	
	
twoway scatter pow_gunmanuf des_gunmanuf_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_gunmanuf des_gunmanuf_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Gun Manuf. (Pow.)") xtitle("Gun Manuf. (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/gunmanuf_si, replace		
	
twoway scatter pow_bigcorps des_bigcorps_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_bigcorps des_bigcorps_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Big Corps. (Pow.)") xtitle("Big Corps. (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/bigcorps_si, replace		
	
twoway scatter pow_bigbanks des_bigbanks_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_bigbanks des_bigbanks_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Big Banks (Pow.)") xtitle("Big Banks (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall))  
	
	graph save Figures/bigbanks_si, replace	

twoway scatter pow_unauthimm des_unauthimm_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_unauthimm des_unauthimm_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Unauth. Immig. (Pow.)") xtitle("Unauth. Immig. (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/unauthimm_si, replace
	
twoway scatter pow_richpeople des_richpeople_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_richpeople des_richpeople_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Rich People (Pow.)") xtitle("Rich People (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall))  
	
	graph save Figures/richpeople_si, replace	

twoway scatter pow_illegalaliens des_illegalaliens_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_illegalaliens des_illegalaliens_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Illegal Immig. (Pow.)") xtitle("Illegal Immig. (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall))  
	
	graph save Figures/illegalaliens_si, replace	
	
twoway scatter pow_police des_police_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_police des_police_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Police (Pow.)") xtitle("Police (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/police_si, replace	
	
twoway scatter pow_military des_military_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_military des_military_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Military (Pow.)") xtitle("Military (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall))  
	
	graph save Figures/military_si, replace	
	
twoway scatter pow_congress des_congress_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_congress des_congress_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Congress (Pow.)") xtitle("Congress (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/congress_si, replace	
	
twoway scatter pow_insuranceco des_insuranceco_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_insuranceco des_insuranceco_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Insurance Cos. (Pow.)") xtitle("Insurance Cos. (Des.)") aspectratio(1) ///
	text(-4 -90 "dependent", size(vsmall)) ///
	text(-4 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/insuranceco_si, replace	


graph combine "Figures/nra_si" "Figures/ice_si" "Figures/whitenat_si" "Figures/gunmanuf_si" "Figures/bigcorps_si" "Figures/bigbanks_si" ///
	"Figures/unauthimm_si" "Figures/richpeople_si" "Figures/illegalaliens_si" "Figures/police_si" "Figures/military_si" "Figures/congress_si" "Figures/insuranceco_si" ///
	,  graphregion(color(white))
	graph save Figures/lowconsensus_si, replace


********************************************
* Figure 6. Distribution of Democrats’ GSDs
********************************************

use "Exported Data Files/sd_party_sdcat.dta", clear

hist des_dem_sd_sdcat2, bin(8) ///
	graphregion(color(white)) xtitle("SD of SD for Deserving (Democrats)") ///
	color(gs10) lcolor(gs5) 

	graph save Figures/hist_des_dem, replace

********************************************
* Figure 7. Distribution of Republicans’ GSDs
********************************************

use "Exported Data Files/sd_party_sdcat.dta", clear

hist des_rep_sd_sdcat2, bin(8) ///
	graphregion(color(white)) xtitle("SD of SD for Deserving (Republicans)") ///
	color(gs10) lcolor(gs5)

	graph save Figures/hist_des_rep, replace

********************************************
* Figure 8. Social Construction Mappings of Immigrant Groups
********************************************

use PartisanWorldviews_Data_Final, clear

	
twoway scatter pow_dreamers des_dreamers_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_dreamers des_dreamers_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Dreamers (Pow.)") xtitle("Dreamers (Des.)") aspectratio(1) ///
	text(-2 -90 "dependent", size(vsmall)) ///
	text(-2 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/dreamers_si, replace	

twoway scatter pow_illegalaliens des_illegalaliens_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_illegalaliens des_illegalaliens_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Illegal Immig. (Pow.)") xtitle("Illegal Immig. (Des.)") aspectratio(1) ///
	text(-2 -90 "dependent", size(vsmall)) ///
	text(-2 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/illegalaliens_si, replace	
	
twoway scatter pow_ice des_ice_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_ice des_ice_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("ICE (Pow.)") xtitle("ICE (Des.)") aspectratio(1) ///
	text(-2 -90 "dependent", size(vsmall)) ///
	text(-2 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/ice_si, replace	

twoway scatter pow_migrantworker des_migrantworker_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_migrantworker des_migrantworker_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Migrant Workers (Pow.)") xtitle("Migrant Workers (Des.)") aspectratio(1) ///
	text(-2 -90 "dependent", size(vsmall)) ///
	text(-2 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/migrantworker_si, replace		
	
twoway scatter pow_refugees des_refugees_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_refugees des_refugees_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Refugees (Pow.)") xtitle("Refugees (Des.)") aspectratio(1) ///
	text(-2 -90 "dependent", size(vsmall)) ///
	text(-2 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/refugees_si, replace	

twoway scatter pow_unauthimm des_unauthimm_si if partisan_rep == 1, ///
	msize(tiny) msymbol(smtriangle) ylabel(none) xlabel(none)  leg(off) mcolor(gs10) ///
	yline(50, lwidth(vthin) lcolor(gs12)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	|| scatter pow_unauthimm des_unauthimm_si if partisan_dem == 1 , ///
	msize(tiny) msymbol(smcircle_hollow) ylabel(none) xlabel(none) leg(off) mcolor(gs2)  ///
	yline(50, lwidth(vthin) lcolor(gs15)) xline(-50, lwidth(vthin) lcolor(gs15)) ///
	graphregion(color(white)) ytitle("Unauth. Immig. (Pow.)") xtitle("Unauth. Immig. (Des.)") aspectratio(1) ///
	text(-2 -90 "dependent", size(vsmall)) ///
	text(-2 -10 "deviant", size(vsmall)) ///
	text(104 -90 "advantaged", size(vsmall)) ///
	text(104 -10 "contender", size(vsmall)) 
	
	graph save Figures/unauthimm_si, replace	

graph combine  "Figures/dreamers_si" "Figures/ice_si" "Figures/illegalaliens_si"  ///
	"Figures/migrantworker_si" "Figures/refugees_si" "Figures/unauthimm_si"  ///
	,  graphregion(color(white))
	graph save Figures/immigration_si, replace


********************************************
* Table 1. Effects of Partisanship on Low Consensus Groups (All Respondents) 
********************************************
use PartisanWorldviews_Data_Final, clear

reg des_nra_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_1
reg des_ice_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_2
reg des_whitenat_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_3		
reg des_gunmanuf_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_4
reg des_bigcorps_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_5		
reg des_bigbanks_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_6		
reg des_unauthimm_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_7			
reg des_richpeople_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep, cluster(fips)
		estimates store tlow_8		
		
* low consensus groups	
	estout tlow_1 tlow_2 tlow_3 tlow_4 tlow_5 tlow_6 tlow_7 tlow_8 using Tables/table_lowconsensus2.txt, replace ///
		cells(b(fmt(3) star label(Coef.)) se(fmt(3) par label(Std. Error))) ///
		stats(N r2 p, fmt(0 2)) ///
		label varlabels(_cons Constant) ///
		starlevels(* 0.10 ** 0.05 *** 0.01) legend ///
		varwidth(25) ///
		modelwidth(12) ///
		mlabels(NRA ICE WhiteNat GunManuf BigCorps BigBanks UnauthImm RichPeople) ///
		collabels(, none) 
		** add F tests manually by looking at the Stata output.		

		
********************************************
* Table 2. Effects of Demographics and Ideological Leanings on Democratic Low Consensus Groups
********************************************			
		
* these are 2 SD larger	
reg des_whitenat_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_dem == 1, cluster(fips)
		estimates store twd_1
reg des_welfarecheats_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_dem == 1, cluster(fips)
		estimates store twd_2	
reg des_sexualharasser_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_dem == 1, cluster(fips)
		estimates store twd_3		
* these are 1.5 SD larger			
reg des_unauthimm_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_dem == 1, cluster(fips)
		estimates store twd_4	
reg des_sexoffender_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_dem == 1, cluster(fips)
		estimates store twd_5		
reg des_illegalaliens_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_dem == 1, cluster(fips)
		estimates store twd_6		
reg des_media_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_dem == 1, cluster(fips)
		estimates store twd_7		


* low consensus groups for Democrats
	estout twd_1 twd_2 twd_3 twd_4 twd_5 twd_6 twd_7 using Tables/table_wedgedem2.txt, replace ///
		cells(b(fmt(3) star label(Coef.)) se(fmt(3) par label(Std. Error))) ///
		stats(N r2 p, fmt(0 2)) ///
		label varlabels(_cons Constant) ///
		starlevels(* 0.10 ** 0.05 *** 0.01) legend ///
		varwidth(25) ///
		modelwidth(12) ///
		mlabels(WhiteNation WelfareCheat SexHarass UnauthImm SexOffender IllegalAliens Media) ///
		collabels(, none) 
		** add F tests manually by looking at the Stata output.		


********************************************
* Table 3. Effects of Demographics and Ideological Leanings on Republican Low Consensus Groups
********************************************		
	
reg des_ice_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 partisan_dem partisan_indep if partisan_rep == 1, cluster(fips)
		estimates store twr_1			
reg des_whitenat_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_2
reg des_nra_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_3		
reg des_bigbanks_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_4		
reg des_richpeople_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_5		
reg des_bigcorps_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_6		
reg des_opioiduser_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_7	
reg des_congress_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_8		
reg des_gunmanuf_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_9		
reg des_insuranceco_dif gender_female race_black race_latino educ_7 age religattend_6  fiscal_conserv5 social_conserv5 ///
	if partisan_rep == 1, cluster(fips)
		estimates store twr_10		
		
* low consensus groups for Republicans
	estout twr_1 twr_2 twr_3 twr_4 twr_5 twr_6 twr_7 twr_8 twr_9 twr_10 using Tables/table_wedgerep2.txt, replace ///
		cells(b(fmt(3) star label(Coef.)) se(fmt(3) par label(Std. Error))) ///
		stats(N r2 p, fmt(0 2)) ///
		label varlabels(_cons Constant) ///
		starlevels(* 0.10 ** 0.05 *** 0.01) legend ///
		varwidth(25) ///
		modelwidth(12) ///
		mlabels(ICE WhiteNat NRA BigBanks RichPeople BigCorps OpioidUser Congress GunManuf InsurCo) ///
		collabels(, none) 
		** add F tests manually by looking at the Stata output.		


		

